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Abstract. Radiant spherical suspensions have an e-periodic distribution in 
a tridimensional incompressible viscous fluid governed by the Stokes-Boussinesq 
system. We perform the homogenization procedure when the radius of the solid 
spheres is of order (the critical size of perforations for the Navier-Stokes 
system) and when the ratio of the fluid/solid conductivities is of order e^, the 
order of the total volume of suspensions. Adapting the methods used in the 
study of small inclusions, we prove that the macroscopic behavior is described 
by a Brinkman-Boussinesq type law and two coupled heat equations, where 
certain capacities of the suspensions and of the radiant sources appear. 
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1 Preliminaries 

One main achievement of homogenization theory was the ability to conceptu- 
ally clarify the relationship between microscopic and macroscopic properties of 
physical systems, at least as far as the periodic approximation could be accept- 
able. The major restriction was the technically impossible interplay between 
different scales: if some quantity varies as the power e" of the size e of the 
mesh, then the case where a < leads to blow up at the limit. This type of 
problems were introduced and solved for the first time by and developed by 
j2j^,4ji^ii^. One major contribution in that direction is the paper by G. Allaire 
13 who clearly underlies the role of critical discriminating scales beyond which 
nothing can be said, but rigidification of elastic systems for instance, and that 
can however generate a transition state where either 'non local' effects [3 E] or 
'coming from nowhere' terms ^ can emerge. 

In this paper, we are insterested in the former case which has been thor- 
oughly explored when non local effects concentrate on rod-like one-dimensional 
submanifolds of the three-dimensional space: see (2] for the Laplacian, for the 
Elasticity system. This geometry enables the formulation of the limit problem 
as a rod-like boundary value problem solved by the density of a Radon mea- 
sure. Our question then was: what happens in other geometries, especially if 
non local effects are to be supported by a cloud of little particles? The physical 
opportunity was the example of thermal flows (see 0121) where highly heat con- 
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ducting spheres are immerged in a Stokes-Boussinesq fluid. It is straightforward 
that for some critical size of the particles (eventually e'^ when the period of the 
distribution is e) the resulting mixture will display a specific behaviour strongly 
discriminating between a trivial case and a classically homogenized case. Our 
concern was then to develop new skills to understand how the expected non 
local effects would be formulated. We found out that the Dirac structure of 
the masses make the classical formulation in terms of a jump term updated and 
that it rather generates an additional source coupled with a capacitary term 
representative of a Brinkman-Boussinesq type law. 

More precisely, the physics of the problem may be described as follows. Solid 
spherical suspensions are e-periodically distributed in a tridimensional bounded 
domain filled with an incompressible fluid governed by the Stokes-Boussinesq 
system. We study the homogenization of the convective movement which is 
generated by highly heterogeneous radiant sources, when the radius of the sus- 
pensions is of £^-order, that is the border case for the Navier-Stokes system (see 
[7]). Assuming that the conductivity and the radiant source of the fluid have £°- 
order, we found that the only regular case in which we have macroscopic effects 
from both the conductivity and the radiation of the suspensions is when they 
are of e^-order. Therefore, we have treated here strictly this case. Nevertheless, 
the present procedure can be easily adapted to the other cases. 

Let n C R'^ be a bounded open set and let 

{k e Z\ Y^^ C ^1} 
The reunion of the suspensions is defined by 

Te ■.= Ukez^B{ek,re), 

where < << e and B{ek, r^) is the ball of radius centered at efc, fc G Z^. 
The fluid domain is given by 

= n\T,. 

Let e*^'^) the last vector of the canonical basis of R'^, n the normal on dfl^ in 
the outward direction and [-Jg the jump across the interface dT^. 

For a > (the so-called Rayleigh number), b > { b (^-^^ denoting the 
ratio of the solid/fluid conductivities), / G Cc{^), g G Cc(ri), where 

Cc(ri) :~ {g G C(il); suppg is compact }, 

we consider the problem corresponding to the non-dimensional Stokes-Boussinesq 
system governing the thermal flow of an e-periodic distribution suspension of 
solid spheres: 
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Set 

14 e iji(rj.;R3), divw = 0}. 

Tlianks to ©, we extend 6^ on by setting 

r = C on T^. 
Tlien, the variational formulation reads: 

V(t;, q) eVeX L^in^), / Vu^ -Vv dx = a / 6'^i;3 da; 

q divu"^ =0 
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fie 

u'^ipV9^ dx = j fipdx + b ( — ) / gipdx. 



, 3 
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We define G H-^{VL) by 



(1) 
(2) 
(3) 
(4) 
(5) 

(6) 

(7) 
(8) 



(9) 



(10) 



gifdx. (11) 



Then, for a > (we shall choose a suitable value for this parameter later), we 
can present the variational formulation of the problem (^-(IHl): 
To find {u^^e^) e 14 X Hl{n) such that 

^{v,^)eV,-KHl{n), (G(u^r),(^;,<^)) =^,(^) (12) 
where the mapping G : x Hl{n) x H^^{fl) is defined by 

{G{u,9),{v,(p)) = a / \7uS/v dx — aa / 9v^ dx 



V6'V(/7 da; + / w(^V6' dx + ^(^— ^ ^ V6'V(/9 dx. 
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In order to prove the existence theorem for problem 112|) . we make use of 
the following result of Gossez. 

Theorem 1.1 Let X be a reflexive Banach space and G : X ^ X' a continuous 
mapping between the corresponding weak topologies. If 

— j — j > oo as \^\x —* oo 

mx 

then G is a surjection. 

Acting as in the proof of Theorem 5.2.2 |H| Ch 1, Sec. 5, we find that the ex- 
istence of the weak solutions of problem H12|l is assured if a is chosen sufficiently 
small. 

Moreover, if is a solution of problem (|12(l . then, by using the weak 

maximum principle, we obtain that 6^ £ L°°{^1), (see Theorem 3.4 |S] Ch 2, 
Sec. 3). 

Remark 1.2 For any a > 0, we have proved the existence of a solution of (|12|) . 
but we do not have a uniqueness result, except if we assume that a > is small 
enough. 

In the sequel, C will denote a suitable positive constant independent of e 
and which may differ from line to line. 



2 Basic inequalities 

Lemma inl and Lemma below are set without proof since it is an adaptation 
of the case p — 2 of Lemma A. 3 and Lemma A. 4 respectively but with 
integrals set on spheres. 

Lemma 2.1 For every < ri < r2, consider: 

C(ri,r2) := {x e R^ ri<\x\<r2}. 
Then, if u G H^{C{ri,r2)), the following estimate holds true: 



iVul 



C(ri,r2) 



> 



wher 



47rrir2 
r2 - ri 



da := 



j- u da — 


j- u da 




1 s., 



47r7" 



• da. 



(13) 



Lemma 2.2 There exists a positive constant C > such that: V(i?, a) G R+ x 
(0,1), Vu e H^{B{0,R)), 



B{Q,B.} 



r>2 

uda\^ dx < C—\Vu\l 
a 



(0,-R)- 
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From now on, wc denote by i?^ a radius with the property << << e, 
that is : 

hm — = hm — = 0. (14) 

Obviously, its existence is insured by the assumption < << e. 
We introduce the measure 

dm^ := — I — I It (x) dx 
47r \rej 

and denote the norm in L'^^ by: 

We denote the domain confined between the spheres of radius a and b by 
C{a,b) := {x e R^ a < < b} 

and correspondingly 

C'=(a,6) :=efc+C(a,6), 
We also use the following notations: 

Ce :=Ufeez,C"-(^e,i?s)- 

S]^^ = dB{£k, R^), Sr^ -.^ Ukez^Sj^^, 
Consider the piecewise constant functions defined after some 6 £ Ho{^) by 

r(a;) = XI ( / Oda]lY>^{x). (16) 



Lemma 2.3 For every 9 E Hq{Q), we have 



[ \e-0^\^ dx 

Jn 


< 




Jn 


I \e-f'\^dx 


< 




[ \vef dx 


[ \e'-f'\^dx 

Jn 


< 


re , 


[ |V6ip dx. 

JCe 



where 0^ and are defined by and U6]) . 
Moreover: 



f \9'\^dx^ f irpdme, / \f'\^dx^ f \f'\^dme 
Jn J Jn J 



(17) 
(18) 
(19) 

(20) 
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Proof. Notice that by definition: 



[ \e-e^\^dx^y] I \e-l eda\'^dx<Y] [ \e-l eda\^dx 



where we have used that 



2 ' 

for every fc e Z^. We use Lemma [2. 21 with 

n = ——, a = 



-V3 
to deduce that 



C— V / |V6'P dx < C— f |V6iP dx 



which shows (|17|) . 

To establish p8|l . we recall the definition: 

( \e-f^\^dx=Y] [ \9-l 9da\^dx 

Applying Lemma 12.21 with R = and a ~ 1, we get the result 

/ le-T^'l^dx < Cr? V /" \Vefdx<Crl[ \Vef dx. 

We come to Indeed, applying Lemma ITTI and (fnji : 

/ V / |/ 9 da -I 9dafdy 

k^^Jy." ^^^^""^ 4^^^^- kti. JC;^^,. 

= C-e^Ir^k^Zll / |v^?|2 da; < C- / \W9\^ dx. 
Finally, a direct computation yields (|20|l . ■ 

Proposition 2.4 For any 6 E iJg(ri), there holds true: 

[ \e\^ dm, <Cmax{l,—) [ \Ve\^ dx. 
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Proof. We have: 



y \0\'^ dm^<2 j \e- f^l^ dme + 2 j |f^|2 dm^ 

= 2 [ \9-f'\^ dm,+2 [ \f^\^ dx 
J Jq 

< Crl [ \Vef dm, +4 [ - dx + 8 [ \e' - dx + S [ \ef dx 
J Jn Jn Jn 

< Crl f-) / |V^|^ dx+C— f |V0|2 dx+C^ f jV^I^ dx+C [ |V^p dx 

<c(— + ^ + l] [ \V9\'' dx<C max (1,—) [ \V9\'' dx 
\re Re J Jn r-e Jq 

m 

Lemma 2.5 For tp € Cc(0) consider the piecewise constant function: 

<fi%x) := V / fdx] lB{ek,r,){x)- 



Then: 



\iin\(p-ip''\m, =0. 

£— >0 



Proof. Notice that 
As we have also 

\B{ek, r,)\ = ^rl card(Z,) ~ H 
then, by the uniform continuity of (p on fi, the result follows. 



3 A priori estimates 

In the sequel, we denote 



r3 

and we assume that 



7. := 5 (21) 



lim7e = 7 e]0,+oo[. (22) 

e— >0 

We denote JTg by 

f Anb f 

T{ip) := fip dx + — gip dx (23) 
Jn <5 Jn 
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Proposition 3.1 We have 

T^^T weakly in H^^{Q) 
Proof. For ip e H^{Q) it follows 

\^ei^)\<\f\nMn.+c(-] \9 







j \9\00 


/ (f dx 







if drrie 



with 



ip dm^ 



1/2 



\ip\'^ dm,!. 



1/2 



Notice that due to (f^ . Proposition 12 . 41 also reads 



I^P dm, < C\Vip\t 



(24) 



(25) 



(26) 



Substituting H25|) and H2()|l into the right-hand side of (|24f) . we get, using Poincare's 
inequality, 

|^,(¥.)| <C|V^|a. (27) 

Now, let (fi € 2?(ri). By the Mean Theorem, there exist e B{ek,r,) such that 

3 



•^e(¥') 



/(/? dx + b 



E 



(7(a;)(y5(x) 
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/ + 5: l!Ir?g(e^)^(f,^) 

/ /^dx + l^^|y/|5(C.1^(e,^). 



There follows 



e I?(r2), limJ>((^)=: / f(pdx-\ / gip dx — J^{ip). 

Jn 3 Jo 

The proof is completed by (|27|l and the density of I?(rj) in 11^(^1). 



(28) 



Proposition 3.2 // {u'^,d^) e x i/o(ri) is a solution of the problem fJii)) . 
flKC? i/ stands for u'^ continued with zero to fl, then we have 



id 9"^ are bounded in Hq[^). 



Moreover, 



6 1 — 1 |V6'^|| < C. 



(29) 
(30) 
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Proof. Substituting v — u'^ in ^ and noticing that 
we get: 

iVw^b, < a|r|o,, (31) 
Seting ip — 9^ in pUI) and taking into account Proposition l3.lL we find 

M'\k+b(^^^ M'\l=J'em<CM^\n (32) 
Noticing that >>1, we deduce from (|^ : 



Therefore 

|Vr|n<C (33) 

and thus 

\0'\n<C. (34) 

Then, follows from Finally, is completed by the estimates (|2U 

and m 

Proposition 3.3 There exist u e H^{Vl;'R?), 9 e Hl{U) and t e L'^{VL) such 
that, on some subsequence, 

^ u in H^{n;R^), 
9^ ^ 9 in H^in), 
^ T in L'^ifl), 
9'^ dm^ ^ T dx in 

where A4f,(ri) is the set of hounded Radon measures on 0, and where denotes 
the weak-star convergence in the measures. 

Proof. From ()29|l . we get, on some subsequence, the following convergences: 

9^ ^9 in Hl{n) (35) 
9' ^9 in L^{n). (36) 
u^^u in iJ(J(r2;R^). (37) 



Moreover, p7|) yields 
which obviously yields 



9'-~e%<c'^^\v9% 



lim|r-r|2 =0. 
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Combining with we infer that 



in L'^iQ). (38) 



We set 



and hence 

0^ dm^ = t'^ dx. 
Taking H29|l and (|26|l into account, we obtain 

dm^ < C. 

We also remark that for any ip £ Cc{^), we have 

ip dm^ —> / ipdx. 



Jn 

Then, using Lemma A-2 of .2,, we find that there exists some r S L^(ri) such 
that, on some subsequence, the following convergence holds: 

dm, ^ T dx, MbiTi). (40) 

Moreover, recall that from (|18|) we have, taking into account (|30|l : 

J |r - f'f dm, < Crl J |Vr|2 dnie < Cr^. (41) 

This implies: 

(r - f") dm, ^ h dx, Mb(n) 
for some h £ and 



< liminf f 16*^ - dm, = 0, 



that is: 

(6*" - f") dm, ^ 0, Al,,(l7). (42) 

Notice that from ifTOIl : 

1^1^ < 2|f^ - + W\l < C-\Ve-\l +C<C, (43) 

and hence, for some f G L^(fi), 

f^^f in L^{n). (44) 
Combining H40() and 142|l . we arrive at 

dm, r dx, A^fc(ri). 

It remains to show that 

T = T. (45) 
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To that aim, let ip S Cc(ri) and let 



We have 



3 / £ 

All \ Tp 



'If dx 



da \lYktfdx 





j- ip dx 




J 9^ip drrip — 






s 









< 


[ {e'' - f^)ip dm. 


+ 


j 




Jn 




Jn 



f^{(p — t/?^) dm. 



From and (|^ . we deduce that 



(46) 



\f'\m=\f'\n<C. 



Moreover, ip G Cc(ri) yields 

\v\m, < C. 

Then, l|4()(l becomes 

I Ut' ~T')pdx\<C\e' ~T'\,r.^+C\p~^'U^. 

Jn 

From 1)41(1 . we infer that 

I / (r^ - f^) p dx\ < Ct, + C\p ~ if'U^. (47) 
Jn 

Thus igZI) and Lemma vield 

lim / (r^ - T')wdx = 0. 

As this holds for every p E Cc{^), the density of Cc{^) in L'^{n) together with 
ij^ and implv that ^ f = r in L^{Q). ■ 



4 The two macroscopic heat equations 

The aim of this section is to pass to the limit as e — > in the variational 
formulation 

V$ei/oi(fi), [ Ve'V'i? dx + b(—\ I VrV$ da;+ 

Jn^ V'"e/ Jt^ (48) 

+ / u'^Vr$ dx = T,{<^>). 
Jn^ 
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Let ip,ip & D(r2) and set 

V^ix) = V f ifi da] lykix), 

rix) - El/ '/^^^ 1 



Let M^"^ denote the fundamental solution of the Laplacian, namely 

AW = in C{r„R^), 
— 1 in r = re, 
W = Q in r = Re. 

The same arguments as in the proof of Lemma A. 3 |^ yield 

re (R. 
{Re -re)\r 



W'{r)^ A — -I] if yeC{re,Re) and \y\ 



Then, we set 

f in ne\Ce, 

w^{x) := I Wix-ek) in Q^ Vfc e 
[ 1 in Te. 

Proposition 4.1 We have 

Proof. Indeed, direct computation shows 
— Ir'' 

('2ir n-K nR., 



£3 V^e ReJ\Re-reJ " (1 - -^) ' 

The proof is completed by (fTljl and if^ . 
For (y5, G 2?(ri), let us define 

$^ = (1- + 

Lemma 4.2 We have 

lim 1$^ - = 0. 

6^0 
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Proof. First notice that ^ in LP'{Vl). Indeed: 

ini Att 

Kln = KkuT. <|CeUT,| = U_i?3 

and linie^^o = by assumption (|14() . As an immediate consequence: 

(I ~ w'')lp -> LP in L^{n). 
Moreover, the uniform continuity of ip over f2 implies that 



so that 



hm - = 



^= = - ^) + ^ in L^{n). 



This achieves the proof. ■ 

Proposition 4.3 // 6*^ is solution of and is given by for any 

ipT^p £ ^{fl), then we have 

lim / • (V$^ + dx 

= j V9 ■ {Vtp + tpu) dx + 47r7 / (6* - t)(V' - ip) dx. 



Proof. First consider 



which reduces to 



/ • (V$^ + ^^m") dx 



Lebesgue's dominated convergence theorem yields V(^lo^\c^ Wip in L^(il). 
Thus, taking (|35|l into account: 

Moreover, 

lln^c."^ - < 1"^ - + l^'Ic.uT, 

and the right-hand side converges to zero because (|37|l yields 

u'^u in L^{n) (58) 

and we apply Lebesgue's dominated convergence theorem to conclude with the 
second term. Thus 

InAc^^^ in L^i^)- (59) 
Now, as e Cc(ri), ((f IsieXCs""^ ^ i^(f^)- Thus, using (|35|l again, 

/ • (plo^\c^u^ dx / V6-ipudx. 
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As a result 

lim 



im / we" ■ {W<i>" + <^>''u") dx^ [ we - {Wip + ipu) dx. (60) 



in,\c. 

Now, we come to the remaining part, namely 



Jce Jce 

Jce 
■■= h+h 



(61) 



We have 



h= I Ve" -{VLp + ifu") dx^ / Vr-V^da;+ / Ve" ■ Lpu" dx. (62) 

Jce Jce Jcs 

In the first term, lc,V<^ -> in L^{n) and Ve" Ve in L^{fl) imply 

• -> 0. (63) 



Ce 

The second term in H62|) is handled by using the estimate: 

\u'\c^ = \lcX\n < \u' -u\n + \u\c,, 

where the right hand side tends to zero due to Using ^ V6' in L^{Q) 
again, we deduce that 

[ we" -ipu" dx^O, (64) 

Jce 

and hence Ii tends to zero. 

It remains to study the integral I2 in (|61|) . To that aim, first notice that 

i2= I we" -Ww^iij" - if) dx = 

Jce 

= (we" -Vw^ii})" ~ip")dx+ fve" -Vw^i^" ~ip)dx (65) 

Jce Jce 

where ip" has been defined by H49|l . The second term in the right-hand side of 
H65|) may be estimated by 

\ [ Ve" ■ Vw"{if ~ if") dx\ < \S/e"\n\Vw"\nW - '/'loo. (66) 

Jce 

As {we) is bounded in H^{Vl), {see Proposition ^Tj, the right hand side of H66|l 

tends to zero by the uniform continuity of Lp over f2. 

Going back to the first term in the right hand side of H65|) , we may write 

Ve" ■ Vw"{il:" - ip") dx 

Ce 

y r d^[\inQdQ r ^ ^r'dr(l ib da ^ I ip 
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reRs 



Y] / i^^\\x-6k\=r, - 0''\\^-ek\=Rj { f ip da - f ipd(7]dai 



{Re - r,) ^ 

from which wc infer that I2 is converging to 

47r7 / (t — 9){ip — if) dx, 

and the proof is completed. ■ 

We are in the position to state a part of our main result: 
Corollary 4.4 The limit (m, 9, r) verifies the following equations: 

uVe~ Ae + 4:nj{e~T) ^ J m (67) 

7(t-6I) = -g in Q.. (68) 
3 

Proof. Consider the variational formulation (|10|l with the test function $ = (f)*^ 
defined by ()57|l for any Lp,%l) £ I?(f7). Then, the left-hand side tends to 



Ve ■ {Vip + ipu) dx + 47r7 / {t - e){ip - ip) dx. (69) 
This is a direct consequence of Proposition ^21 together with the remark that 

VrV$'' dx = 



since is constant on every B{ek,r^), k E Z^. 

The convergence of the right-hand side is obtained by using the uniform conti- 
nuity of ^ and by Proposition Thus we find the variational formulation of 
and the proof is completed. ■ 



5 The homogenized problem 

Proposition l3.3l vields the existence of some u G H^{n; 'Rp) with div(u) = and 
for which the following convergence holds on some subsequence 

u^-^u in H^{fl;R^). 

From [7], we find that there exists an extension of the pressure (denoted by p^) 
and some p G L'^{fl) such that 

p^^p in L^{n)/Ti. 

We denote by (w^,q^) G H^{C{r,, §)) x Ll{C(r^, §)) the only solution of the 
following Stokes problem 

divw,*^ = in C{r,,^), 
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-Aw^ + Vq^ = in C(r„|), 
= if r = r^, 

= e^"') if r = -. 
^ 2 

Consequently, we define 

r if xeTe, 

i e« if xen,\U,ez^C'{re4)- 

For g we set u = (pv^ in (jSl and then using the energy method like in 

^ we find the equation that the velocity field satisfies in H^^{n): 

- Au + Gn^u = -Vp + aOe'-^'^ in Q. (70) 

Finally, we summarize the results of Proposition 13.31 CoroUarv 14. 41 together 
with the relation H7U|I into our main theorem. 

Theorem 5.1 If{u^,p^) is a solution of problem then the following con- 

vergences hold on some subsequence 

9' ^ 9 in H^in), 
9^ dm^ T dx in A4b{^), 

where (u,9) G i7Q(fl;R'^) x HqIQ), which stand for the macroscopic velocity 
and temperature of the fluid, and r S L'^{VL), which stands for the macroscopic 
temperature of the vanished suspensions, form a solution of the following system: 

divM = in f2, 

—An + 67r7u = — Vp + a9e''^^ in fi, 
uVe ~ A9 + ATT-f{9 -t) = f in Vl, 
Anb 

An^ij — 9) = g in Vl. 

3 

Remark 5.2 In the present case, with suspensions of critical size, the Brinkman- 
Boussinesq equation was an expected result; nevertheless, our proof is different 
from that of [Jj, which treated the homogenization of the Navier-Stokes equa- 
tions for perforated domains in a similar case. 

Remark 5.3 Our two-temperature model, with 7 as transfer coefficient, is the 
macroscopic effect of the assumption on the the ratio of the fluid/solid conduc- 
tivities. 
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Remark 5.4 The appearance of the source term —^9 in the second macro- 
scopic heat equation is strictly the consequence of the assumption on the mi- 
croscopic radiation. 
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